Statistical Model Predictions for Particle Ratios at y^s/vJv = 5.5 TeV 
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Particle production in central Pb-Pb collisions at LHC is discussed in the context of the Statistical 
Model. Predictions of various particle ratios are presented with the corresponding choice of model 
parameters made according to the systematics extracted from heavy-ion collisions at lower energies. 
The sensitivity of several ratios on the temperature and the baryon chemical potential is studied in 
detail, and some of them, which are particularly appropriate to determine the chemical freeze-out 
point experimentally, are indicated. We show that the p/p ratio is most suitable to determine the 
baryon chemical potential while the f2~/K~ and Q~ /t7~ ratios are best to determine the temperature 
at chemical freeze-out. 
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I. INTRODUCTION 

Particle production observed in heavy-ion collisions al- 
lows a systematic study of the thermal properties of the 
final state. In a wide energy range, from the SIS up to 
RHIC, the yields of produced particles have been shown 
to be consistent with the assumption that hadrons origi- 
nate from a thermal source with a given temperature and 
a given baryon density. These properties have been quan- 
tified by comparing the measured particle ratios with 
Statistical Model calculations. By using only two ther- 
mal parameters, a successful description of particle ratios 
measured in heavy-ion collisions over a wide range of cen- 
ter of mass energies could be made ^. The extracted 
chemical freeze-out parameters, the temperature T and 
the baryon chemical potential fis can be characterized 
by a constant average energy per hadron {E)/{N) of ap- 
proximately 1 GeV 13 • The extrapolation of this freeze- 
out curve towards vanishing ^b is bounded by the criti- 
cal phase transition temperature as calculated in Lattice 
Gauge Theory j^]. 

In view of the success of the Statistical Model, we dis- 
cuss in this paper the expectations for particle ratios in 
central heavy-ion collisions at the LHC energy. The cal- 
culations are performed at the chemical decoupling point 
extrapolated from the freeze-out systematics. The parti- 
cle ratios are discussed with respect to their sensitivity to 
variations in the thermal parameters. Based on this anal- 
ysis, we suggest observables which are best suited to ex- 
tract experimentally the freeze-out conditions in heavy- 
ion collisions at the LHC energy. 

At LHC energies, due to the large number of pro- 
duced particles, a novel type of data analysis will be 
possible. In particular, the thermal origin and chemi- 
cal decoupling conditions in heavy-ion collisions can be 
studied for various event classes and even within a sin- 
gle event. Consequently, the freeze-out conditions can 
be extracted event-wise, leading to distributions of the 



temperature and baryon chemical potential determined 
event-wise. Instead of a full analysis within the Statisti- 
cal Model, the proposed particle ratios can be utilized to 
extract these values. 

Furthermore, at these energies most particles might be 
produced via hard collisions which could lead to particle 
ratios substantially different from those of the Statistical 
Model. Deviations from these predictions might lead to 
new insight into the hadronisation mechanism. 



II. STATISTICAL MODEL 

The Statistical Model and its application to heavy- 
ion collisions has been summarized recently [l|. In the 
present calculations we apply the Statistical Model with 
completely equilibrated flavor production. Due to the 
expected large multiplicity of strange particles in cen- 
tral Pb-Pb collisions at LHC energies, we use the grand- 
canonical formulation of strangeness conservation. Thus, 
strangeness, baryon number and electric charge are de- 
termined by corresponding chemical potentials. It should 
be noted that for very peripheral heavy-ion collisions, the 
suppression of the strange-particle phase space due to 
canonical effects might occur even at LHC energies. 
The charge chemical potential is constrained by the ini- 
tial isospin asymmetry of the Pb nuclei, whereas the 
strange chemical potential /ig, depending on both T and 
fiB, is determined by strangeness neutrality. Thus, any 
particle ratio is uniquely determined by only two param- 
eters, T and /is at chemical freeze-out. 

In order to estimate particle ratios in Pb-Pb collisions 
at LHC, one needs to determine the expected range of 
thermal parameters. For this purpose we apply the pa- 
rameterized Y^s^vTV-dependence of fis as it was recently 
extracted from the Statistical Model analysis of hadron 
multiplicities in heavy-ion collisions [J|. The baryon 
chemical potential fiB, extrapolated to the LHC energy. 



TABLE I; Particle ratios in central Pb-Pb collisions at freeze- 
out conditions expected at the LHC: T = (170±5) MeV and 
Hb = ill MeV. The given errors correspond to the varia- 
tion in the thermal parameters. Additional, systematic un- 
certainies in the ratios of the right column arise from unknown 
decay modes. They are smaller than 1% in general, but reach 
3% in the E~ /A ratio and 7% in the p/vr" and the A/p ratios. 
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is, according to this parameterization. 



MB(\/sJviv = 5.5 TeV) 



1 MeV. 



(1) 



Similar estimates of ur at the LHC energy have been 
obtained recently in 0) la S 13 • 

Considering the variation of the chemical freeze-out tem- 
perature with coUision energy in heavy-ion coUisions, ev- 
idently T is an increasing function of ^/snn- However, 
from SPS to RHIC the change of the chemical freeze-out 
temperature is only moderate. Within the statistical and 
systematic errors the temperature T « 170 MeV is con- 
sistent with the value derived at the top SPS and RHIC 
energies [H H ES [H El . Lattice QCD results on the 
critical temperature, extrapolated to the chiral limit and 
at vanishing baryon density, concur within statistical er- 
ror with this value |3(|. Thus, we use this value as an 
estimate of the chemical freeze-out temperature at LHC. 
However, to account for possible uncertainties in the ex- 
trapolations of the chemical freeze-out parameters to the 
LHC energy, we consider a possible shift of T by ±5 MeV 
and a variation of fiB between and 5 MeV. 

From the sensitivity of different particle ratios on 
freeze-out parameters we deduce the reliability of the 
model predictions and the ability to determine the pa- 
rameters experimentally. To identify the relevant parti- 
cle ratios that can be used to constrain experimentally 
the chemical decoupling conditions in heavy-ion collisions 
at the LHC energy, we consider a broader range of pa- 
rameters. The ratios are studied for 150 MeV < T < 
180 MeV and for < /xb < 30 MeV. 



III. PARTICLE RATIOS AT LHC 

With the values of the thermal parameters expected at 
LHC and discussed in the previous Section we calculate 
a set of particle ratios using the Thermus package [13 ■ 
The predicted particle ratios for Pb-Pb collisions at LHC 
are summarized in Table HI and ITll The errors arise from 



TABLE 11: Left column: Resonance/stable particle ratios in 
central Pb-Pb collisions at freeze-out conditions expected at 
the LHC: T = (170±5) MeV and fin = itt MeV. The given 
errors correspond to the variation in the thermal parameters. 
Additional, systematic uncertainies arise from unknown decay 
modes. They are smaller than 1% in general, but reach 2% 
in the <^/K" ratio and 3.5% in the P°(770)/tv~ ratio. Right 
column: Resonance width, taken from 
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the variation of the temperature and of the baryon chemi- 
cal potential, as indicated. Additional systematic uncer- 
tainties, not included in the Tables, are due to incom- 
plete knowledge of decay modes of resonances into sta- 
ble particles. We estimate that due to the uncertainty 
in branching ratios, the p/tt" and the A/p ratios might 
vary by up to 7%, whereas for the S~/A ratio a change 
of 3% is possible. The other ratios are less sensitive to 
the resonance properties and affected by less than 1%. 
Within the considered range of thermal parameters, the 
antiparticle/particle ratios are insensitive to any varia- 
tion in contributions from resonances, due to the decay 
chains being equivalent for particles and antiparticles. 

One of the expected features of particle production at 
LHC is that there should be a rather negligible differ- 
ence between the yield of particles and their antiparti- 
cles. This is a direct consequence of the low net-baryon 
density expected in the collision fireball at LHC at chem- 
ical decoupling. With the range of the chemical potential 
< /J,B < 5 MeV used in the actual calculations, the ra- 
tios of particle to antiparticle yields, called h/h hereafter, 
are indeed seen in Table Q] to be close to unity. Figure Q] 
illustrates the dependence of h/h ratios on T (left) and 
on hb (right). The calculations at fiB ~ 27 MeV serve 
as cross checks with experimental data. At RHIC, at 
y'sNN = 200 GeV, several particle ratios were measured 
in Au-Au reactions by the BRAHMS 14], PHENIX [3, 
PHOBOS O and STAR |17| collaborations. Good 
agreement is observed for all ratios under study. As an 
example the p/p ratio is indicated in Figure ^ 

On the other hand, the ratios of particles with differ- 
ent masses and different quantum numbers are predomi- 
nantly controlled by the freeze-out temperature and the 
particle masses (see Sec. IIVI below). The values of such 
yield ratios vary strongly with the particle mass differ- 
ence, as seen in the right panel of Table Q] 
Some of the ratios are expected to change only slightly at 
LHC from their values measured at the SPS and RHIC. 
This can be seen from Fig. 12 showing little sensitivity 
of these ratios with T and fiB ■ A particular example are 
the K/tt ratios which, within errors, are consistent with 
those measured at lower energies. Large deviations from 



SPS and RHIC measurements are seen for particle ra- 
tios built with baryons and antibaryons. This is a direct 
consequence of a substantial decrease in baryon stopping 
at LHC, which, in the Statistical Model, results in a de- 
creasing baryon chemical potential. 

Particle densities in central Pb-Pb collisions at 
y/J/VTV = 5.5 TeV have been recently calculated within 
relativistic hydrodynamics |1^. Their results on ratios 
of pions, kaons and protons are comparable to our find- 
ings given in Table ^ However, the predictions of this 
hydrodynamic model depend strongly on the value of the 
decoupling temperature. A similar sensitivity on the de- 
coupling temperature is found in both models. 

Resonance yields were suggested as a sensitive probe 
of the fireball expansion dynamics [13 ■ The resonances 
might be sensitive to collective effects and interactions 
with the surrounding medium. Such interactions can 
result in resonance broadening or mass shifts [20I l2ll| . 
Assuming that resonances are produced according to a 
thermal distribution, their initial multiplicity can be cal- 
culated. Deviations from the thermal yields can be at- 
tributed to processes which take place in a later stage 
of the evolution after chemical freeze-out. In this late 
stage, the expansion dynamics might be driven mostly 
by elastic and pseudo-elastic hadron scattering. Also 
the reconstruction of short-lived resonances which de- 
cay in this stage will fail if at least one daughter par- 
ticle undergoes an interaction with the medium. On 
the other hand, pseudo-elastic processes might regenerate 
resonances |23| and consequently, the yield of resonances 
at thermal freeze-out can differ from the one calculated 
in the Statistical Model at chemical freeze-out. This has 
been observed in measurements of resonance yields in 
heavy- ion collisions at RHIC [22| ■ The observed yields of 
short-lived resonances like p or K* were found to differ 
from the model predictions even by a factor of two [l|. 
Benchmarks are essential to possibly evaluate modifica- 
tions in resonance yields. Therefore, resonance/stable 
particle ratios, as expected at chemical freeze-out, are 
calculated in the Statistical Model without including dy- 
namical effects in the final state. They are displayed 
in Table ^ together with the width of the resonances. 
Yields of shorter lived resonances, i. e. those with a large 
width, are expected to be affected more strongly [2^ . 



IV. SENSITIVITY TO FREEZE-OUT 
PARAMETERS 



In this section we discuss the dependence of different 
particle ratios on the values of the thermal parameters. 
The main objective of this study is to identify observables 
that could serve as sensitive experimental probes for the 
chemical freeze-out conditions at the LHC energy. In this 
context we consider the antiparticle to particle ratios as 
well as the ratios of strange and non-strange particles of 
different masses and quantum numbers. 
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FIG. 1: Antiparticle/particle ratios i? as a function of /is for 
T = 170 MeV (left) and as a function of T for /is = 1 MeV 
(right). The horizontal line at 1 is meant to guide the eye. 
The p/p ratio (averaged over the data of the 4 RHIC exper- 
iments at y/SNN ~ 200 GeV) is displayed (gray horizontal 
line) together with its statistical error (gray band). As illus- 
trated, hb ~ 27 MeV (dashed line) can be read off the Figure 
directly within the given accuracy (vertical gray band). 



A. Antiparticle to Particle Ratios 

Figure ^ shows the dependence of different antiparti- 
cle/particle ratios, h/h, on /zb (left) and T (right). For 
vanishing baryon chemical potential, the density of par- 
ticles is identical to that of antiparticles, thus h/h — 1. 
For finite and increasing baryon chemical potential, the 
h/h ratios are decreasing functions of fis- Such proper- 
ties of antiparticle/particle ratios are qualitatively well 
understood in the Statistical Model. The mass terms 
in the single particle and antiparticle partition functions 
are identical. Consequently, the leading dependence of 
the h/h ratios on fiB is determined by. 



h/h (X exp[-2(S^B + S'ms)/^], 



(2) 



where B and S are the baryon and strangeness quantum 
numbers of the particle, respectively. In the expression 
above, feed-down contributions from resonance decays 
are ignored but they have been included in the model 
calculations. However, due to equivalent contributions of 
resonance decays to particles and antiparticles, Eq. ^ 
provides a good description of h/h ratios at the LHC 
energy. In the antibaryon/baryon ratios, the /iB-term 
dominates over the /xg-term because strangeness conser- 
vation implies that the baryon chemical potential is, in a 
strangeness neutral system and at a fixed temperature. 
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FIG. 2: Ratios R of particles with unequal strangeness con- 
tent as a function of /is for T — 170 MeV (left) and as a 
function of T for hb = I MeV (right). 



always larger than the strangeness chemical potential. 
This, together with the opposite sign of the quantum 
numbers B and S of strange baryons, results in a weaker 
sensitivity on ^.b of h/h ratios for hadrons with increas- 
ing strangeness content. Consequently, as seen in Fig. ^ 
at a fixed value of ^b , the antibaryon/baryon ratios are 
increasing with strangeness content of the hadrons. 

The K~/K+ ratio shows qualitatively a similar de- 
pendence on the baryon chemical potential as the an- 
tibaryon/baryon ratios. This is due to the strangeness 
content of the kaons and the fact that fis is an increasing 
function of /i^. The K~/K+ ratio has no explicit depen- 
dence on the baryon chemical potential and is therefore 
always larger than the antibaryon/baryon ratios at the 
same fiB and T. 

For small values of fj,B , as expected in heavy-ion colli- 
sions at the LHC energy, the h/h ratios are only weakly 
dependent on temperature, as seen in Fig. ^ (right). 
Thus, these ratios do not constrain the chemical freeze- 
out temperature in heavy-ion collisions. On the other 
hand, the high sensitivity of h/h on fiB makes such ratios 
ideal observables to quantify /is experimentally. As seen 
in Fig. 1, this is particularly the case for the p/p ratio 
which exhibits the strongest dependence on the baryon 
chemical potential. 



B. Mixed Particle Ratios 

From the preceding Section one can conclude that ra- 
tios of particles with equal masses are not good ther- 



mometers of the medium. To extract the chemical freeze- 
out temperature in heavy-ion collisions, one should con- 
sider ratios that are composed of hadrons having different 
masses. The dependencies of such ratios on /jb and T are 
illustrated in Fig. [3 Within the range of the freeze-out 
parameters considered in this figure, there is a rather 
weak sensitivity of the ratios on T and fiB ■ In the case of 
the rj^/S^ ratio this can be seen analytically from the 
expression. 
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where, in this simplified formula, feeding from resonance 
decays is ignored, with m^- and m^- denoting the cor- 
responding particle masses. 

There is no explicit dependence on the baryon chemical 
potential in all ratios considered in Fig. |21 as also seen in 
Eq.lPJ. Changes in /is affect the ratios only indirectly 
through the /is-dependence of the strangeness chemical 
potential. The values of fis are nearly proportional to /is 
in the range studied here. However, for all values of ^b 
considered in Fig. 2, it is always smaller than 8 MeV. 
Consequently, the dependence on the baryon chemical 
potential is comparatively weak. 

The temperature dependence of the fJ^/S^ ratio is dom- 
inated by the mass term, since the difference ((toq- — 
m-^-) ~ 350 MeV) is much larger than the strangeness 
chemical potential (/ig < 1 MeV for ^b ~ 1 MeV and for 
all considered values of T). However, to correctly quan- 
tify the temperature dependence of the ratios of particles 
with different masses, one needs to include contributions 
from resonance decays. 

Figure |3| (upper and middle panel) displays particle 
densities from different sources. Particles called primary 
are directly produced from a thermal system with given 
temperature and chemical potentials. The second contri- 
bution, called decay in Fig.jSj originates from short-living 
resonances which feed the yield of stable hadrons. The 
sum of both contributions results in the final particle den- 
sity. 

There is no feeding to 17" hyperons from heavier reso- 
nances. In contrast, only about 50% of the S~ hyperons 
are primary, while the other half originates from reso- 
nance decays. Since both contributions exhibit a similar 
temperature dependence, as demonstrated in Fig. 13 the 
sensitivity of the Vt~ /'B.~ ratio on T can still be estimated 
from Eq. (PJ. The ri~/S~ ratio increases by about 1% 
per 1 MeV temperature increase. 

For light hadrons like kaons and pions, the feed-down 
contributions exceed the thermally produced yield. As 
seen in Fig. |21 (lower panel) for the thermal conditions 
at LHC, approximately 75% of the pions are expected 
to be produced by resonance decays. Furthermore, the 
contributions from resonance decays feature a tempera- 
ture dependence deviating from that of primarily created 
hadrons, as seen in Fig. 13 the former exhibit a steeper in- 
crease. Therefore, the temperature-driven increase in the 
heavy/light hadron ratios is progressively diluted when 
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FIG. 3: Upper and middle panel: Particle densities of 7r~, 
H~, K~ and Jl" hadrons as a function of T for fiB ~ 1 MeV. 
Different sources are indicated. The thermally produced, pri- 
mary hadrons together with those from decays sum up to the 
finally observed densities. The primary and final density of 
Q~ hyperons is identical, (i.e. there is no contribution from 
resonances) . 

Lower panel: Fractional amount of tt", K~ and H~ hadrons 
from primary production and from decays of the final densi- 
ties, shown as a function of T for hb ~ 1 MeV. 



hadrons with lower masses are considered in the denom- 
inator. In the case of the K/tt ratios, the contribution to 
the total pion yield from decays compensates fully for an 
increase in the kaon yield with increasing temperature. 
Consequently, the K/tt ratios are almost independent of 
temperature. 

This observation allows a stringent test of the model. If 
the experimental values for the K/tt ratios deviate signifi- 
cantly from the predictions, they can hardly be reconciled 
with the Statistical Model. This would be an indication 
for a non-equilibrated system at chemical freeze-out. 

From the discussions above and in the previous Sec- 
tion, it becomes clear that neither antiparticle/particle 
ratios nor ratios composed of hadrons with small mass 
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FIG. 4: Particle ratios R involving hyperons as a function of 
T for /is = 1 MeV. 



differences are well suited to deduce the chemical freeze- 
out temperature. However, ratios with hyperons, in par- 
ticular the hyperon/pion ratios, are excellent observables 
to extract the chemical decoupling temperature in heavy- 
ion collisions. Figure ^ shows the temperature depen- 
dence of some particle ratios containing hyperons. The 
Statistical Model predictions of these ratios are shown up 
to T = 200 MeV. 

In general, an upper limit for the chemical freeze-out tem- 
perature at LHC is given by the value of the QCD critical 
temperature at vanishing chemical potential. A recent 
lattice QCD finding i2^ indicates that in (2-|-l)-flavor 
QCD the critical temperature can be even as large as 
195 MeV. This value is significantly larger than previ- 
ously expected ^. Figure ^ illustrates that the particle 
ratios are sensitive to temperature variations in the range 
that includes the above value of the critical temperature. 
However, at temperatures larger than 190 MeV, the Sta- 
tistical Model results start to suffer from an imprecise 
knowledge of high-mass resonances and their decay prop- 
erties. 

The sensitivity of different particle ratios on temperature 
is quantified in Fig. 0] by the variation of the ratio R per 
1 MeV change in T relative to the value of R calculated 
at T = 170 MeV. The results given in this figure show 
that the temperature dependence is strongest for parti- 
cle ratios with a large mass difference. This statement is 
also valid for ratios composed of antiparticles. 
The largest (2.5%/MeV) sensitivity on temperature is 
exhibited by the fl~ /t:~ ratio. However, the variation is 
significantly smaller than that expected from only ther- 



mally produced particles. This is because at LHC only 
about 25% of all pions are expected to be directly pro- 
duced, while the remaining fraction originates from bary- 
onic and, predominantly, from mesonic resonances [25J- 
The feed-down contribution to pions increases with tem- 
perature and reduces the sensitivity. 
In view of the large decay contributions to the pion yield, 
the ri^/K~ ratio might be a better thermometer of the 
medium, since the feeding to kaons is noticeably smaller. 
Figure 0] demonstrates that the two ratios, ri~/7r~ and 
f2~/K~ , are similarly sensitive to temperature variations. 
The combination of the larger kaon mass and smaller 
feeding compared to pions results in comparable temper- 
ature dependencies. 

The 'Er l'K~ ratio has a similar mass difference as the 
f2~/K~ ratio. However, the very large feeding of heavy 
resonances to pions destroys the expected similarity in 
the temperature dependence of these ratios. The sen- 
sitivity of the hyperon to proton ratio on temperature 
is rather weak due to the small mass difference of con- 
stituent particles. 

From Fig. ^it is obvious that the 17~/7r^ and fi^/K^ 
ratios are best suited to determine the chemical freeze- 
out temperature. Experimentally, the latter might be 
more easily accessible and more precisely known, since 
the contributions from weakly decaying particles, which 
have to be obtained with high accuracy, are larger for tt^ 
than for K^ mesons. 



V. CONCLUSION AND SUMMARY 

Predictions of the Statistical Model for different parti- 
cle ratios for Pb-Pb collisions at the LHC energy are pre- 



sented. The sensitivities of various ratios with respect 
to the temperature and the baryon chemical potential 
as well as the contribution of resonances are discussed 
and analyzed. We have shown that the p/p ratio is the 
best suited observable to extract the value of the baryon 
chemical potential at chemical freeze-out. The 17^ /tt^ 
and the ri~/K~ ratio are proposed as thermometers to 
extract experimentally the chemical freeze-out tempera- 
ture in central Pb-Pb collisions at LHC. 



Hard processes could dominate the particle production 
at LHC energies. From these interactions, very different 
particle ratios than those of a thermal source might be 
expected to be created. Deviations from the predictions 
given here will help to understand the dynamics and, 
possibly, the time scale of the hadronisation procedure. 
Furthermore, in events selected e.g. by high-energy jets, 
the resulting particle ratios might be different, indicat- 
ing perhaps a higher freeze-out temperature or perhaps 
deviating completely from the thermal picture. Such se- 
lection criteria might make it possible to differentiate be- 
tween early and late freeze-out. On the other hand, the 
frequent interactions between the constituents close to 
the phase boundary Sq could lead to a complete reshuf- 
fling and a statistical distribution might be observed at 
the end, independent of the preceding history. 
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